Backscattering: an overlooked effect of General Relativity? 
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The total flux of outgoing radiation in a strong gravitational fleld decreases due to backscattering 
if the sources are close to an apparent horizon. It can cause detectable changes in the shape of 
signals. Backscattering could well be of relevance to astrophysics and would constitute a new test 
of the validity of general relativity. An explicit bound for this effect is derived for scalar flelds. 
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Backscattering can be described as a phenomenon 
which causes the radiation of a massless field to disperse 
outside those null cones defined by the initial impulse, ft 
is known in older mathematical literature as the 'break- 
down of Huyghens principle' |^ and it is interesting for a 
number of reasons. First, as a non-Minkowski spacetime 
effect, it offers a new way of testing general relativity. 
Second, it can be important in astrophysics, both in or- 
der to explain the radiation coming to us from regions 
adjoining a black hole and to infer information about the 
sources of the radiation. Third, a strict upper bound on 
the magnitude of the backscattering effect, such as we 
derive here, offers numerical relativists an independent 
test of the correctness of their numerical codes. 

Our model is a spherically symmetric massless scalar 
field emitted from a region near to, but outside, a 
non-rotating black hole. These simplifying assumptions 
should not seriously restrict the validity of our conclu- 
sions. For example, the propagation of electromagnetic 
fields, we believe, should obey inequalities similar to the 
ones derived below. These estimates, while they break 
down close to the horizon, allow us to distinguish the re- 
gion in which the backscattering may play a significant 
role from that in which it is of no importance. They offer, 
to our knowledge, the first quantitative measure of this 
strong field effect. Others who work on backscattering 
adopt quite different approaches 

We consider a foliation of the spacetime by using the 
polar gauge slicing condition, trK = K^; that is, with a 
diagonal line element 

ds^ = -l3iR,t)j{R,t)dt^ + ^^^dR^ + R^{r,t)dn^ , 

j[R,t) 

(1) 

where t is the time, i? is a radial coordinate which co- 
incides with the areal radius and dil^ — d9^ + sin^ Odip'^ 
is the line element on the unit sphere, < < 27r and 
< 6 < tt; P and 7 go to -1-1 at infinity. 

For a massless scalar field the stress-energy is T^,y — 



V^0V ^(f) ~ gfiijW a4>^°'(j>/'2. The matter energy density is 
p — — Tq and the matter current density is J = —Tqr/ (3. 
(do ± jOr) are the outgoing and ingoing null directions. 
We define radiation amplitudes 

h+{R,t) = h{R,t) = ^{-do+jdR){Rcj)) (2) 
27 

h.{R,t) ^ h{~R,t)^ ^(do+jdR)iR^). (3) 
27 

One can show that 

;3(i?) = e-«^/ri^[('^+-^)^+(''-^)1 , (4) 

and 

^(i?) = 1 _ ^ + ^ _ f3ir)]dr , (5) 

where h — Jj^ dr[h-^-{r) + h^{r)] = ^0 . 7 can be 
expressed in the following useful form |^] 

7W - (1 - ^ + m (6) 

where toq is the asymptotic (total) mass and uiext is a 
contribution to the asymptotic mass coming from the 
exterior of a sphere of a radius R, 
1*00 

m,xt{R) = 47r 1 [[h+{r) - h]^ + [h^{r) - h]^) dr . 

(7) 

The scalar field equation W^d^cj) = can be written as 
a single first order equation on a 'symmetrized' domain 
-c» < i? < c» [§ by writing h{R) = h+{R) and h{-R) = 
h^{R) as 

{do+jdR)h^{h-h)^^ . (8) 

Eq. (^, together with the definitions of h, h, (3, and 7, is 
equivalent to the Einstein equations coupled to the scalar 



1 



field. The external mass changes along an outgoing null 
cone according to 



[do + ldR)me^t{R) = -Sir-r^h^ - hf 



(9) 



The polar gauge allows us to express the metric directly 
in terms of the matter, Eqs. (^- 1|) , where all the integrals 
are in the exterior region. The local and global Cauchy 
problems for the above system are solvable in an external 
region bounded from the interior by a null cone ( [H, [^). 

Following Eq.(|), {-f-P)/R ^ -2m{R)RP/\RfwheTe 
■m{R) = mo — mextiR) is the Hawking mass at a ra- 
dius R. Eq.(||) gives a 'red-shift' due to the h term on 
the right-hand-side (determined by the mass function, 
m{R), rather than the Schwarzschild mass, mo) and a 
'backscattering' due to the h term. 

Let us define 
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2m(i?)' 
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2m{r)dr 



r^{l — 2m{r)/r) 



(10) 



where the integral is taken along an outgoing null ray. 
This allows us to rewrite Eq.(0) as 



(ao + 75fl)(l 



2m{R) 
R 



)h 



- 2m{R) 2m{R)R(3 



R 



(11) 



Thus (1 — 2fri{R)/R) is the redshift factor and the right- 
hand-side of Eq. ( prj ) determines the backscattering. 
It is natural to write Eq.(H) as 



and by using Eq.(^) we get that 

/? , 1 ^ 1 
77" r — 2m{r) ~ r — 2mo 



hydr 



(12) 



(13) 



The factor P/'yr can be taken out of the integral in 
Eq.(p^ and replaced with its value at r = i?. The re- 
mainder is then essentially right hand side of Eq.(^. Let 

us choose an e and an Ra such that mextiRA)/mo < e 
and Ra > 2mo(l + e). For any R > Ra we obtain 



1 > l3iR) > > e 



0{mext/emo) . 

(14) 



In the same vein, using Eq.(j6|) one gets j{R) ~ 1 — 
2mo/i?-l-0(mea;t/emo). Thus the effect of the matter on 
the geometry can be controlled. Notice that in regions 
sufficiently close to the horizon, even if mej-^/mo ^ 1, 
a small cloud of matter can still strongly influence the 
geometry. In what follows, however, we will restrict our 
attention to the region outside R ~ 3mo, since only in 
that region of spacetime can we get sensible analytic es- 
timates. This is equivalent to choosing e > 1/2. 



For R > 3mo the scalar wave equation Eq. (11) can be 
approximated as 

rfi , « \f^ 2mo - 2mo 2moi? 
{do+jdR){l —)h^h{l — ) , (15) 



R 



R ' \R\^ 



with the error terms of order rrif.xt/'^o- In the limit of 
'mext/iTT'O <C 1 the equation describes a scalar field prop- 
agating on a fixed Schwarzschild background. 

One can show, improving a coefficient in an inequality 
of @, that 



\h\< 



y/mext 



i?l/2 /8^(1 _ 2mo ) 



(16) 



(A similar estimate with mo rather than rUext, appears 
in Q.) Let us introduce the Regge- Wheeler coordinate 
§ r* = i? + 2mo ln(i?/2mo - 1) so that -/dn = dr- . The 
solution of Eq. ( [l5| ) can be estimated above and below by 
solutions of 



(5„ + a,o((i-^)/.) = 



Eqns.(n7^) are solved by 



2mo(l - i^)7^ 
i?5/2^87r(l- ^) 



(17) 



l-^Ah{T\t) = h^{r* -t) 



ir,t) (1-^) 
(r-t,0) i?,Vfl:-2m( 



■dv 



(18) 



here the dv in Eq.(|l8|) is dr* , with dR/dr* — I — 2mo/R. 
If we ignore the second term then the standard redshift 
is obtained. 



hir*{0)+T, t) 



/ 1 _ _2mo_\ 
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2mo 



/i(r*(0),0), (19) 



fl(r*(0)+T)^ 

while the integral in Eq.(|l8|) can be solved to give 

dR. 



IGTrmo 



J^^mo f ^fl = 

V 27r UJ Ri^{R-2mo 

^ (1 - ^) + arctan 



R.-2mo 



2mo 



(20) 



Thus the total backscattering is bounded by a term of 
order y^m^It/mo- In the limit where ■mf.xt/rriQ is small, 
we recover the usual gravitational redshift along an out- 
going null ray with an error of order ^Jriiext/'mQ- This 
holds even in a strong gravitational field and no 'quasi- 
static' assumption need be made. We need not assume 
R S> mo, but only R > 3mo. 

Let us assume that initial data at t = represents a 
pure outgoing wave. In other words, h-(t = 0) = 0. 
We select an Ra > 3mo on the initial slice such that 
"niextiRA) / fno ^ 1- Finally we assume, in addition to 
^ at infinity, that = for R < Ra- This guaran- 
tees that the radiation is bounded away from the black 
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hole. Consider the future outgoing lightray from {Ra , 0) . 
This is well outside any event horizon, which would be at 
approximately R = 2mo, so this lightray is really outgo- 
ing all the way to null infinity. Note that mext{RA,0) 
is the maximum value of ruext over the whole wedge 
bounded by the initial slice and the outgoing null ray. 
We will estimate /i_ and h along this null ray. The inte- 
gration of Eq. (^) along this ray yields an estimate of the 
total energy flux across this surface in the inward direc- 
tion. This will be the total energy loss from the outgoing 
wave due to backscatter. 

Choose a point on the null ray from Ra and label it 
by (i?i,Ti). To calculate at this point, consider the 
ingoing future null ray which passes through this point 
and integrate Eq.(|^ along this ingoing ray. This will 
start from the initial hypersurface at some point (i?2,0) 
with R2 > Ra- Along this null ray R monotonically 
decreases while rriext monotonically increases. 

To get an explicit estimate, the integral in Eq.Jl^) 
can be further approximated; since R monotonically de- 
creases along the ingoing lightray, we can replace the 
y/1 - 2TOo/i? by y/1 ~ 2mo/Ri. This yields 



Hi 



dR 



R2 flVH-2mo 



< 



-Ri r dR ^ 



3R, 



Ri 



3/2 Y (_R,i-2mo) 



Thus we arrive at 



fnext 



3 V 2TTm, 



3/2 

Ry/R - 2mo 



(21) 



(22) 



We can write, using Eq.(||) and (p2|). 



m. 



3/2 



i?(i?-2mo)i/2 



(23) 



and a similar inequality with a minus sign to give a lower 
bound of h. 

These equations are integrated along the outgoing null 
ray from Ra to give 



\Rh{R,t)\ < \rHRa,0)\ + IJ'^mo 



dR 



[R - 2mo)3/2 
(24) 



Since we demand that 0(i?/i,O) ~ 0, the first term in 
Eq.(p4r) vanishes. Thus we can bound h by 



4 / rrte^t 



3/2 

3 \l 2irmo i?.(i?.A-2mo)i/2 
3/2 



3 Y 2ii-mo ii(fl-2mo)i/2 



(25) 



The last term in Eq.(p5|) is strictly larger than as 
given by Eq.(p2|) if i? > Ra > 4mo (i. e., when all 
radiation is placed outside 4too). Thus 



\h-h-\ < \h\ + \h_\ < 



A I 3/2 r> 

4 / ruext Ra 



3 V 27rmo R{Ra - 2mof/^ 



(26) 

Therefore, from Eq.(|^) the total change in nig^xt satis- 



fies 



16/2mo\ ( l~mo/RA , , , 



From this expression we can see how sensitive the 
amount of backscattering is to the location of the in- 
nermost null cone. This estimate becomes meaningless 
if i?2 ~ 3.5mo because /S.mext > Weat- However if 
Ra = 6r7io we get that less than 25% of the total en- 
ergy in the exterior field is backscattered. In the case of 
a neutron star, where 2mo/R < 0.1 (on the surface of 
the star), we have an upper bound for the backscattered 
energy of 2% of rriext ■ 

In deriving the bound Eq. ( p7|) , a number of truncations 
and approximations were used. While most of them are 
sharp in the sense that configurations exist that turn the 
inequality into an equality, we do not believe that all of 
them can be sharp simultaneously. Therefore Eq.(p7|) is 
clearly an overestimate. We present it here, not because 
it is the best that can be done with this technique, but 
because it is simple to derive, it is in analytic form, and 
is easy to interpret physically. 

For example, in Eq. ( pl| ) we ignored the term depending 
on i?2. We can include this term if the initially outgoing 
pulses are far enough from the apparent horizon [Q, so 
as to sharpen the estimate to niext 



^ ^ W /2mo\2(l-mo/i?2) 



(28) 



where ~ 0.4. Applying this to the case of neutron 
stars, when the conditions assumed above hold true, we 
find that the maximal amount of backscattered radiation 
cannot exceed 1 percent. 

Eq.(p^ can be integrated in closed form, (Eq.(po|)), 
but in Eq.(pT|) we approximated it to get a simpler ex- 
pression because we needed to integrate the result again. 
Numerical integration would obviate the need for this ap- 
proximation entirely, but the result would be much less 
transparent. 

The tail term can also be bounded using this method. 
The initially outgoing field, h^, generates a weaker ingo- 
ing field, which enters the 'no-radiation' zone behind 
the wavefront. This, in turn, scatters again off the grav- 
itational field to generate a new outgoing field, which 
turns up at null infinity at a later time, after the first 
burst of outgoing radiation has gone. This is the so- 
called 'tail' term. Eq.(^5|) gives us an estimate for h in 
the 'no-radiation' zone. Substituting this back into the 
scalar wave equation (^|) , an estimate of the second or- 
der along the outgoing null ray is 
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\-fh+{R,t)\ < 



4 / TOe^t ™0 

3 y 27rmo VflA-2mo J-R-A 



\/J?.A-2mo 



dR 



In the limit as i? ^ cxd along the null cone we get 



|/i+(oo,oo)| < - 



27rmo ^1 - 2mo/i?2 



(29) 



(30) 



This should be compared with the leading term in (18), 
/iQ. Using the definition of the external mass mext, we 
get that the tail term is smaller than the leading term by 
a factor {rrio/Ro)^- 

We expect that photons will behave in a way similar 
to the massless spin zero field that we have discussed 
here. Gravitational physics predicts two phenomena for 
a radiating plasma surrounding a black hole or a neutron 
star: Redshift diminishes the intensity and frequency of 
the outgoing radiation while the total energy in the ra- 
diation (as measured by the mass function) remains un- 
changed. Backscattering, the other effect, weakens the 
overall outgoing radiation, so that the total energy that 
reaches infinity is reduced. In addition, backscattering 
changes the shape of extended signals - the leading part 
weakens while the rest of the impulse gains in intensity. 
This conclusion follows from a careful analysis of (||) and 
is supported by the numerical results of |Q, which show 
that up to 10% of the radiation emitted at i? = 3mo is 
'shifted' inside the main impulse; even if it reaches infin- 
ity, it does so with a significant delay. 

There are two astrophysical situations where those two 
different aspects of backscattering will play a role. The 
net efficiency of the black hole - matter system is reduced 
below what is expected when only the redshift factors are 
taken into account. This may be of significance in mod- 
elling 'too faint' galactic nuclei fuelled by black holes, 
such as the nuclei of M87 ||]. Second, backscattering may 
leave imprints in X ray bursts resulting from energetic 
processes on the surface of a neutron star or in collisions 
of compact bodies, such as neutron star - neutron star, 
neutron star - black hole, or black hole - planetoid merg- 
ers. Backscattering will deform the peak and contribute 
a tail term to the radiation emitted from such short-lived 
sources. Gamma ray bursts (GRB) are believed by some 
pof to arise from such collisions and should reveal traces 
of backscattering. The absence of these effects would give 
great support to the fireball scenario of GRB's 0. 

Another more immediate application is in numerical 
relativity. Much work has been done in constructing 
codes to analyse the Einstein - massless scalar field model 
[ p2[ . Bounds on the magnitude of the backscattering 
effect, such as derived here, would offer numerical rel- 
ativists a reliable test of the correctness and long-time 
stability of their codes. 
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